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ABSTRACT 


This thesis investigates the fundamentals of multiple—input—multiple—output 
(MIMO) radio communication systems with space-time codes. A MIMO system was de- 
signed using the Alamouti space-time code. The modulation technique was binary 
phase-shift keying (BPSK). Matlab with Simulink was used to simulate the design, 
which was tested in both an additive white Gaussian noise (AWGN) channel and in a 
multipath fading channel with AWGN. Theoretical performance was derived for both 
channels and compared to simulated results. The original receiver design was changed to 
incorporate a maximal-ratio combiner (MRC) receiving technique with channel state in- 
formation (CSI). The theoretical performance for this design was determined and com- 


pared to simulated and published results. 
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EXECUTIVE SUMMARY 


In modern radio communications, the demand for high speed, reliable communi- 
cation together with efficient use of the spectrum and power are the prime technical crite- 
ria for communications engineers. However, these implementations must operate in a se- 
vere environment. They must confront rich scattering effects and multipath propagation 
under limited power constraints. This not only restricts the speed of communications but 
the reliability as well. Multiple—input—multiple—output (MIMO) communication systems 
with space-time coding provide increased reliability and increased capacity without 
bandwidth expansion. MIMO communications systems have application to existing and 
future military radio communications systems where channel bandwidth is restricted but 
increased capacity is necessary. If the information rate needs to double and must fit 
within the existing bandwidth, then MIMO communications systems can provide the in- 


crease in capacity. 


This thesis investigated the fundamentals of MIMO systems with space-time cod- 
ing. A simulation of a simple design was created to demonstrate and analyze perform- 
ance. The space-time code chosen was the Alamouti scheme using binary phase—shift 
keying (BPSK) with a two transmit and two receive antenna MIMO system. The system 
design was implemented in Matlab with Simulink. The design was tested in an additive 
white Gaussian noise (AWGN) channel and then a multipath fading channel with 
AWGN. The original design was modified using a maximal-ratio combiner (MRC) tech- 
nique and tested in a multipath channel with AWGN. The simulated performance results 


are compared with theoretical analysis as well as published results. 


The results showed that the designed MIMO system performed within expectation 
of theoretical analysis in AWGN. It outperformed a single-channel BPSK system, need- 
ing 3.0 dB less power for equal performance. However, the designed system failed in a 
multipath channel. The MRC design was incorporated in the receiver. As a result, the 
new design was able to communicate successfully over the multipath channel. The per- 
formance of the new design matched expected theoretical performance as well as pub- 


lished results. The designed MIMO system with the MRC receiver was able to achieve 


XV 


full diversity order with the Alamouti space-time code which was consistent with pub- 
lished results. It was able to achieve a significant advantage over a single-channel BPSK, 


or single—input—single—output (SISO) system, in a multipath channel with AWGN. 


XV1 


I. INTRODUCTION 


A. BACKGROUND 

In modern radio communications the demand for high speed, reliable communica- 
tion together with efficient use of the spectrum and power are the prime technical criteria 
for communications engineers. For example, the implementation of wireless Internet 
with video conferencing or other media applications requires wideband, high data rate 
and reliable communication systems. However, these implementations must operate in a 
severe environment. They must confront rich scattering effects and multipath propaga- 
tion under limited power constraints. This not only restricts the speed of communica- 
tions, but the reliability as well. Added to this is the limited available spectrum that is al- 
located for wireless systems. Limited bandwidth means limited data rates, or speed of 
communications. What if a system can be implemented that provides reliability and in- 
creased capacity within limited bandwidth constraints? This would give the communica- 
tions engineer more avenues for trade-offs in design. The demand for increased capacity 
within limited bandwidths is certainly relevant to existing and future military communi- 
cations systems. A multiple—input—multiple—output (MIMO) communication system 
with space-time coding holds such a promise for engineers. “Recent research in informa- 
tion theory has shown that large gains in capacity of communication over wireless chan- 
nels are feasible in multiple—input—multiple—output systems.”[1] A review of MIMO re- 
lated literature finds the quote or variations of it several times. The primary focus of this 


research was to investigate the fundamentals of MIMO systems. 


B. GOALS AND METHODOLOGY 

The goal of this thesis was to examine MIMO systems with space-time coding 
and their relevance to military communication systems. In order to do so, a simulation of 
a simple system was created to demonstrate and analyze the design’s performance. The 
space-time code chosen here is the Alamouti scheme using binary phase-shift keying 
(BPSK) with a two transmit antenna and two receive antenna MIMO system [1]. The 
modulator and simple receiver design were created and implemented in Matlab with 


Simulink. The design was realized in equivalent baseband form to facilitate comparison 
1 


with published theory. The design was tested first in additive white Gaussian noise 

(AWGN). Theoretical performance was determined and compared with simulated results 
of a single channel BPSK system, or single—input—single—output (SISO) system. The de- 
sign was then progressed to a multipath fading channel with AWGN. The simulated per- 


formance results were compared with theoretical results as well as published results. 


C. BENEFITS OF STUDY 

Space-time coding in MIMO systems is new technology. There is a great deal of 
research effort in this technology [1]. The demand for high speed reliable communica- 
tions systems in a limited frequency spectrum for commercial use is mimicked by mili- 
tary applications. However, some military applications also need capacity and conceal- 
ment. If increased capacity and performance at a certain power level is achieved, then 
the power level may be reduced by sacrificing capacity. Nevertheless, this may be at a 
capacity that is already acceptable. The reduced power level may lead into a region of 
low probability of detection communications and concealment. Notwithstanding, MIMO 
and space-time coding holds promise for additional trade-offs for the communications 


engineer. 


D. THESIS ORGANIZATION 

This thesis is organized in four remaining chapters. Chapter II introduces the 
fundamentals of the MIMO model. It discusses capacity, the channel and the significant 
aspects of MIMO and space-time coding. Chapter III introduces the design of the 
MIMO communications system using the Alamouti space-time code [1]. The transmit- 
ter, receiver and logic of the receiver’s creation will be discussed. Examination of the 
theoretical and simulated performance of the MIMO system design is detailed in Chapter 
IV. It progresses through performance in AWGN to performance in a multipath fading 
channel with AWGN. Chapter V reviews the results and concludes with recommenda- 
tions for future study. Additionally, there are two appendices: Appendix A is a list the 


acronyms, and Appendix B contains the block parameters used in the simulations. 


YW. MULTIPLE-INPUT—MULTIPLE—OUTPUT 


This chapter lays the groundwork for understanding and analyzing MIMO sys- 
tems in a multipath fading environment. The MIMO model is introduced and the MIMO 
capacity equation is derived. The capacity of a MIMO system is greater than that of a 
SISO system [2]. The increase in capacity of MIMO systems is the driving force of 
much of the research in this area. Next, an overview of the multipath fading channel is 
presented. The signal will be subjected to multipath fading, so a baseband model is de- 
scribed to facilitate analysis. Then, diversity techniques and the maximal-ratio combiner 
(MRC) are described. A firm understanding is necessary for analysis of the performance 
of the design. Finally, space-time codes are introduced. Space-time coding is a mecha- 
nism to approach the capacity of a MIMO system [1]. The analysis in this thesis is per- 
formed in the baseband transmission model. Thus, a definition of the baseband transmis- 
sion model is introduced first. This chapter provides the fundamentals needed to analyze 


the design of the MIMO system with the Alamouti scheme. 


A. BASEBAND TRANSMISSION MODEL 
The baseband transmission model is the lowpass equivalent of the bandpass 
transmission model [3]. It consists of the complex lowpass representation or baseband 


equivalents of bandpass signals and the bandpass channel transfer function as illustrated 


in Figure 1. In Figure 1, the input bandpass signal is represented by x(t) and its Fou- 
rier transform X ( ri ) . The bandpass system’s impulse response is represented by h(t) 
and its Fourier transform H(f). Lastly, the output of the system is represented by y(t) 
and its Fourier transform Y ( f ) . Let the baseband equivalent of each of the signals and 
transfer function be distinguished by a tilde, 1.e., X(t) ; h (t) and y(t) . Similarly, the 
Fourier transforms of the baseband equivalent system are represented by X(f), H(/) 


and aga 





Band—Pass 
t @) 
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x) = ¥(f) 
3 Complex = 
*(t) Low—Pass — 
£(s) ea PCs) 











Figure 1. | Bandpass and equivalent lowpass system (After Reference [3].). 


The input signal x(t) is related to its baseband equivalent by [3] 
x (i) = Rel x (tje (2.1) 
where Re[ | represents the real part of the expression in the brackets and /, is the car- 


rier frequency of the bandpass system. The baseband equivalent signal X(t) can be rep- 


resented by its inphase and quadrature components x,(¢) and x, (t) by [3] 


X(t) =x, (t)+ jxo(t), (2.2) 
where 
x,(0)=Re[¥(1)], 2.3) 
and 
Xp (t) = Im[ ¥(¢) |. (2.4) 


The Im| | is the imaginary part of the function inside the brackets. 


The magnitude of the frequency response of the bandpass system H ( f ) and that 
of the complex lowpass equivalent system H ( Fa ) are illustrated in Figure 2. The base- 


band equivalent transfer function H ( ra ) can be obtained by shifting the positive fre- 


quency components of the bandpass transfer function ( £ ) by f, to be centered about 


zero [3]. The transfer function is assumed to flat over the bandwidth B. 
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Figure 2. Bandpass and equivalent lowpass transfer function (After Reference [3].). 


From Figure 1, the Fourier transform of the output signal of the bandpass system 


Y(f) and baseband equivalent system is given by [3] 
Y(f)=X(Sf)A(S). 23) 


The baseband equivalent is given by 


Y(f)=X(f)A(S). (2.6) 
It is assumed that the bandwidth of X(t) is small enough such that the magnitude and 
phase of H ( f ) can be considered constant in this bandwidth of interest; specifically, 

H(f)=he”, (2.7) 
where /h is a positive real number and @ is in the range of [0, 27]. Then h (t) is given 
by 


h(t) =he!’d(t). (2.8) 


where 5(t) is the Dirac delta function. The baseband equivalent output signal j)(t) is 
the convolution of X(t) and h(t), given by 


ioe) 


H(t)=] 8 


t) h(t —t)dt 
) 


he!’ d(t—t)dt (2.9) 


The bandpass output signal y(t) is related to the its complex equivalent j(t) by [3] 
y(t)=Re| j(t)e?*" J. (2.10) 


Like the input signal x(t), the baseband equivalent output signal »(t) can be decom- 


posed into inphase and quadrature components by [3] 


H()=y,(t)+ io(t), (2.11) 
where 
y,(t)=Re[ 5(0)], (2.12) 
and 
¥o(t)=Im[¥(4)]. (2.13) 


The bandpass system can be obtained from its baseband equivalent. The informa- 


tion of contained in x(t) and y(t) are preserved in ¥(¢) and (t) [3]. In the rest of this 


thesis, the signals and systems are represented by their baseband equivalent. 


B. THE MIMO MODEL 
The MIMO model is described in Figure 3. The development of the model dis- 


cussed in this section is largely based on Reference [1]. The model consists of N, trans- 


mit and N, receive antennas. From each transmit antenna, a signal x,, where 


j=1,2,3,...,N,, is transmitted. The signal x, represents the baseband equivalent ¥,(t). 


(The function of ¢ has been dropped for convenience.) Each signal x, passes through 
N, channels represented by /,,. It is assumed that the bandwidth of the signals x, is 
small enough that Equation (2.9) holds true. The symbol /, is a complex number given 
by 


ef”, (2.14) 








and will be referred to as the channel coefficients. At each receive antenna, white Gaus- 


sian noise (AWGN) n,, where 7 =1,2,3,..., N,., is added to the signal. The symbol n, 
represents the AWGN baseband equivalent 7, (t) [1]. The received signal is represented 


as r., where i=1,2,3,...,N., and is the baseband equivalent of the received signal 7, (t) ; 


Transmitter Receiver 











Figure 3. | MIMO model (After Reference [1].). 


At any receive antenna /, the received signal 7, is the sum of the transmitted sig- 


nals through the channel and the AWGN and is given by [1] 


= MyX thx, +... +hiy Xy +N; 


(2.15) 


Thus, the system as a whole can be represented as a system of linear equations given by 


[1] 








r hy, hy, LN, x, n, 
7 h h,,--*h x n 
2 21 22 2N, 2 2 
mile . a roe ae (2.16) 
™y, | Aw 1 hy hy N, Xw, Ny, 
or, in matrix form, as 
r=Hx+n (2.17) 


where r represents an NV, x1 column matrix, H is an N, x N, matrix, x isan N, x1 col- 
umn matrix, and n isan N, x1 column matrix [1]. 

Let each element in x be a zero-mean, independent identically distributed (IID) 
Gaussian random variable [1]. From information theory, the optimum distribution for 


transmitted signals is Gaussian, thus they are assumed to be Gaussian random variables 


[1]. The autocorrelation matrix R.. of the transmitted signal x is defined by [1,4] 


R SE xx, (2.18) 


xX 


where A”’ is the complex conjugate transpose of the matrix A. Since x is a vector of 
zero-mean IID Gaussian random variables, they are uncorrelated [4,5]. The correlation 


matrix is diagonal. The covariance of x is given by [1,4] 
C.. =R,,. (2.19) 


The total average power P of vector x is given by [1] 


p=e| Sis; ]-Dellsf ]-m(R.) (2.20) 


8 


where tr(A) is the trace of matrix A, defined as the sum of its diagonal elements. It is 


assumed that each transmit antenna transmits equal power. The power from any one an- 


tenna is given by [1] 


2 P 
P=E| lx FS (2.21) 
Then R.. is 
P 
R,.. Sap tt . (2.22) 


t 


The symbol I,, represents the identity matrix of size NV . 


It was assumed in the previous section that the transmitted signal’s bandwidth x, 


is narrow enough that the frequency response of the channel is flat. Also we assume that 
there is no attenuation due to path loss or amplification due to antenna gain except as 
measured in the channel matrix H. Furthermore, we assume that the channel representa- 
tion matrix H is fixed during a symbol period, i.e., each element of H is a fixed coeffi- 
cient. In the discussion of a multipath fading channel to follow, the appropriateness of 


these assumptions will be made clear. 


Let the normalization constraints be [1] 


N, 
Ly 
j=l 


°=N, fori=1,2,...N,. (2.23) 





The sum of the magnitudes of each channel coefficient leading into any one receive an- 


tennais N,. Ifthe channel coefficients are modeled as random variables, then the nor- 


malization constraint will apply to the expected value of the left hand side of Equation 


(2.23) [1]. This normalization is necessary for comparison to a SISO system. 


The matrix H constitutes the channel state information (CSI). It is assumed that 
there is perfect knowledge of the CSI at the receiver. While this assumption is never true, 
it can be approached. For example, each coefficient can be estimated by sending training 


signals to the receiver [1]. 


The AWGN at each receive antenna is described by a single column matrix n. 
The elements of n are zero-mean, equal variance, IID complex Gaussian random vari- 


ables [1]. As with the transmitted signals, the autocorrelation matrix n is given by [1,4] 


R,, =E[nn™ |. (2.24) 


nn 


Similar to the transmit matrix x , since each element of the noise is assumed zero-mean, 
the covariance matrix is equal to the autocorrelation matrix. Each element of n is an IID 
random variable with equal variance; thus, they are independent and uncorrelated [5]. 
The power of each noise element is equal to the variance. The autocorrelation function is 


[1,4] 
R,,, = Oly, . (2.25) 


Let the average power at each receive antenna in the absence of noise be repre- 


sented by P. which is given by [1,4] 


Doo (2.26) 











where the channel coefficients have been assumed fixed, i.e., not random variables. The 
( y represents the complex conjugate. Substitution of Equations (2.21) and (2.23) into 


Equation (2.26) gives 


E| [x P [=n a (2.27) 


R= >It 
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Then the average signal-to-noise ratio (SNR) at each receive antenna, represented by y, 
is given by 
P 
VaR (2.28) 
0, 


It was assumed that there was no path loss or gain through the channel with the exception 


of that measured by the channel coefficients. The total received power P,, becomes 


tot 





N,. 
P= >) P=N,P=tr(R,,). (2.29) 


The received signal autocorrelation matrix is 
R, = E[rr™ | : (2.30) 


Using Equation (2.17) in Equation (2.30), we obtain 


' (2.31) 
= b| (Hx+n) (Hx) : +n”) | 
Expansion of Equation (2.31) gives 
R,, = E| Hx(Hx)” +Hxn"” +n(Hx)” +n” | 
(2.32) 


=E| Hx(Hx)" |+E[Hxn” |+B| (Hx)"n]+E[ nn” | 


since the expectation of a sum of random variables is the sum of expectations. The trans- 


position of a product of matrices has the identity [6] 
(Hx) =x7H". (2.33) 


Then Equation (2.32) can be written as 





R,, = E[ Hx" H” |+5] Hxn ]+ 5] x Hn |+E/ nn” J. (2.34) 
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If H is deterministic and x and n are independent of each other, then they are also un- 


correlated; thus, 











E| Hxn"” |= HE[x]E| n°” | 
= * * * a * (2.35) 
E| x""H'n|=E| x” |H"E[n]. 
Since both random variables x and n have zero-means, Equations (2.35) becomes 
E| Hxn” |=0 
oetee (2.36) 
E| x"7H'"n | =0, 
and Equation (2.34) reduces to 
R,, = HE| xx” JH” +E/ nn” J. (2.37) 
Substitution of Equations (2.18) and (2.24) into Equation (2.37) gives 
R,,=HR, H’+R,,. (2.38) 


Finally, there are two equations that define the MIMO model, Equations (2.17) 
and (2.38): 


r= Hx+n, 


: 2.39 
R, =HR,H’+R,. Cy 


The first is in terms of random vectors and the second is in terms of the autocorrelation 


matrices. With this information, the capacity of a MIMO system can be derived. 


Cc. CAPACITY OF MIMO SYSTEMS 
The derivation of the capacity of MIMO systems presented here is roughly drawn 


from that found in Reference [1]. First, consider Equation (2.17) 
r=Hx+n. (2.40) 
Any Nx WN, matrix H can be decomposed by the singular value decomposition (SVD) 
theorem into a product of matrices [4] 
H=UDV”. (2.41) 
12 


For reference, the matrix Q is defined by [1] 


0 HH” when NV, < N, (2.42) 
H’H when N, < N.. | 


Note that the singular values of H are defined to be the square roots of the eigenvalues of 


Q [1,4]. Let the eigenvalues of Q be represented by 4, where i=1,2,3.... Now, Equa- 


tion (2.40) can be rewritten as 
r=UDV’x+n. (2.43) 


The matrices D, U and V have the following properties [1,4]. Matrix D isan Nx WN, 


matrix of nonnegative, real, singular values of H. Let the singular values of H_ be repre- 


sented by p,. Then the matrix D is given by 








p, 0 3 0 
0 P2 
D=/0 0 Pn, when N, <N,, (2.44) 
0 0 0 
| 0 0 0 | 
Py 0 
0 Pp, aus 0 aa4 
D=|. a ; e when NV, < N,, (2.45) 
0 0 Py, 0 
P 0 0 
0 p> 0 
D=}. when NV,=N,=N. (2.46) 
0 0 Pn 
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Making the assumption that 


Pre Py ice pees (2.47) 


it is seen that some of the higher numbered p, can be zero. If the rank of H is r, then 
there are 7 nonzero singularities [4]. The matrix U is an N,xN, unitary matrix com- 


prised of the left singular vectors of H [1,4] 
U=|u, u, be uy |- (2.48) 


The left singular vectors of U , represented as u,,u,...u, , are defined as the eigenvec- 


tors of HH” [1,4]. The matrix V is an N,N, unitary matrix comprised of right singu- 


lar vectors of H given by 
V=lv, Vv, ue Vie |s (2.49) 


The right singular vectors of V , represented as v,,v,...v, are defined as the eigenvec- 


tors of HH [1,4]. Lastly, a unitary matrix A has the following properties [6]: 
(2.50) 


and 
AA™ =AA™ =I, (2.51) 
Both matrices U and V have these properties. 


Let the following transformations take place on r, x and n [1]: 


x SV ox; (2.52) 
U 


Then 


r= (U" y r’, 
x=(V7) x, (2.53) 
n= (U" y n’ 
Substitution of Equations (2.53) into Equation (2.43) gives 
(U7) r'=UDV" (V7) x'+(U7) a’, (2.54) 


Multiplying both sides of this equation by U’’ and using the properties defined in Equa- 
tion (2.50) and (2.51), we obtain 
U7 (U7) =UTUDV" (V7) x4 U7 (U7) 
Ir’ = IDIx' + In’ (2.55) 
r’=Dx'+n’. 


Therefore, Equation (2.55) is an equivalent representation of Equation (2.40). 


Expanding Equation (2.55) when N, < N,., we get 


























iH vA 0 ° ny 
a 0 Vr, “ 0 : ny, 
: i : x, 
Rodeo £0 Fe lice eae am (2.56) 
Ty, +1 0 0 pti : My My. 
I J 2 8 mf Lak | 
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Similarly, when NV. < N,, Equation (2.55) expands to 


‘ JA, 0 wets 
rn |_|0 ae i 0 


N, 0 0 


Transmitter 








Figure 4. 
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_ oS 














(2.57) 





The equivalent MIMO model N, < N,. (After Reference [1].). 


Transmitter 





Figure 5. The equivalent MIMO model N. < N, (After Reference [1].). 


Note that the 0 in Figures 4 and 5 means that no signal passes in that channel. For either 


situation, V,< N, or NV, <N,, the equivalent received signal at any receive antenna, 


where r is the rank of H, is 
r= JAxital where i=1,2,...r. (2.58) 


The rank of H can be no larger than the minimum value between N, and N. 
[1,4,6]. Furthermore, from Figures 4 and 5 and Equation (2.58), it can be seen that the 
equivalent received signal at any receive antenna 7’ is a direct couple between the 
equivalent transmit signal x/ with an applied gain factor VA; , a singular value of H. No 


longer are the received signals a linear sum of all transmit signals multiplied by a channel 


coefficient, as was the case in the general MIMO model. 


The equivalent channel describes r direct uncoupled SISO channels. The capac- 
ity of the system as a whole is just the sum of the individual capacities [1]. Shannon’s 


capacity equation is [7] 


C=Whlog,[1+SNR], (2.59) 
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where C is the capacity in information bits/second, W is the bandwidth of the channel in 
Hertz and SNR is the signal-to-noise ratio. Thus, the capacity of the MIMO system is 
given by 


C=WY log, [1+SNR/], (2.60) 
i=l 
where SNR! is the signal-to-noise ratio for the i-th equivalent channel. To determine the 
maximum spectral efficiency C/W , it is sufficient to find the equivalent received signal 


and noise power. The equivalent received power in the i-th channel is the power of the 
transmitted signal multiplied by the associated eigenvalue. In order to obtain the power 


of x’, the autocorrelation matrix is found; that is, 
R,, =E| xx"? |. (2.61) 
The substitution of Equation (2.52) into Equation (2.61) gives 
R,.=E [vx (v"x)" | (2.62) 


Using the transposition identity defined in Equation (2.33), we see that Equation (2.62) 


reduces to 


R,., =E| V7xx''V | 
=V7E| xx"? |V (2.63) 
=V"R.V. 


Substitution of Equation (2.22) in Equation (2.63) gives 


R.=V" 1, Vv 
Nt 


xx 


(2.64) 
Z —VTL,V. 


t 


The matrix V is an N,N, unitary matrix; using the identity in Equation (2.51), we see 


that Equation (2.64) reduces to 
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P 
R. ara (2.65) 


Similarly, to find the power of the equivalent noise the autocorrelation of n’ is 


used, which is defined as 
R,, =E[n'n”” |. (2.66) 


Once again, substituting the transformation defined in Equation (2.52), we get 


“E 


R,,,=E [un (U"n) | 
=E/ Unn'U | (2.67) 
=U"R,U. 


The substitution of Equation (2.25) into Equation (2.67) gives 
R,, =U"o,1,,U. (2.68) 


The matrix U isan N, x WN, unitary matrix; using the identity in Equation (2.51), we see 


that Equation (2.68) reduces to 


n 


Roy = Oly, - (2.69) 


Therefore, the SNR of the i-th equivalent channel is 








SNR! = ane : (2.70) 
195 
Using Equation (2.70) and (2.60), we find the capacity of the MIMO system to be 
A.P 
C=WY) log,|1+— 
2, &> Noe | 
(2.71) 


=W log, +25 


i=l t-o 
where J, are the eigenvalues of Q and r is the rank of H, assumed to be given by 


r=min(V,,N,). 12) 
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The function min(V,, V,) represents the minimum between NV, and N,. 


Ideally, Equation (2.71) should be in terms of the channel matrix H. Let the ei- 
genvalue equation be defined by [1,6] 


(AL, -Q)y =0 where y #0. (2.73) 
The matrix Q is defined in Equation (2.42). 


0 HH” when N, < N, (2.74) 
H’H when N, < N.. , 


Equation (2.73) has a non-trivial solution if 4 is the root of the characteristic equation 
[6]. The characteristic equation is 
p(A)=det(AI, -Q)=0. (2.75) 
The eigenvalues of a square matrix Q are the roots of the characteristic equation 


p(A); therefore, the solution to p(A) can be written as [6] 


P(A)=(-A)A-A)-- (4-4, =0 


r 2: 
=| [U-4)=0, : ie 


where 4,,/,,...,4, are the eigenvalues of Q. Combining Equation (2.75) and (2.76), we 


get 


[[(@-4) =aet(21, -Q). (2.77) 


i=l 
This is true for all values of 4 in the real domain [1,4]. 


Noting the similarities of the left hand side of Equation (2.77) to (2.71) and mak- 
ing the substitution defined by [1] 


Azote (2.78) 
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in Equation (2.77), we get 


. No? No? 
es 9 et a te 
He = a! oad, Q 


i=l 











(2.79) 
- A,P P 
1+——) = det} I, + . 
at No? [ N, x) 
Therefore, the equation for the capacity of a MIMO system becomes 
z AP 
C=Wlo 1+— 
&> ni Noe 
(2.80) 





P 
C=W log, [ae + No 0}} 


It is appropriate to discuss three examples to illustrate the result. For the first ex- 
ample, let the system be a SISO system, i.e., a system with a single transmit and a single 
receive antenna. From Equation (2.59), the normalized capacity, i.e., the maximum spec- 


tral efficiency, becomes 


< =log, [1+ SNR]. (2.81) 


Let the SNR be 10.0 dB. Thus the normalized capacity is 3.46 bit/s/Hz. 


For the same example but this time Equation (2.80) is used. In this case 


Q=1,I, =1 and N, =1. Using this information in Equation (2.80), we get 


C=W log, a + 4) 
Oo, 


C P 
ry +4) (2.82) 
& 


=log,(1+SNR). 
W £2 ( ) 
Thus, Equation (2.80) reduces to Equation (2.81). 


For the next example, consider a MIMO system with the same number of transmit 


and receive antennas, that is, N, = NV. =N. Let the channel matrix H_ be defined by [1] 
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[Ne OO esc TT 
H=|° VN Y _ : (2.83) 
OG. -@. OD a AR] 








This is equivalent to N independent, uncoupled SISO channels. The matrix Q is de- 


fined by HH and given by, 


N 0 0 0 
_|0 N 0 0 
ee m (2.84) 
0 0 0 N 
= NI,. 


Using Equation (2.80), we find the normalized capacity to be 


C P 
— =log,| det) I, +—— . 2.85 
Substituting Equation (2.84) into Equation (2.85), we get 
C NP 
W = log, fax, ea) 
‘Pe 
-ve,{de(1, +251, } (2.86) 
0, 
= log, faa{(- 3] I }} 
0, 


Since the determinant of a diagonal matrix is the product of the diagonals [8], Equation 


(2.86) becomes 


N 
C P 
— =log, t + 5 
W 2 
ae (2.87) 
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The capacity of the channel using a MIMO system as defined by matrix H in 
Equation (2.83) gives a capacity N times that of the SISO case. This makes sense since 
the channel is equivalent to N independent SISO channels. The capacity increases line- 
arly with the number of antennas. If N =2, then the normalized capacity for SNR of 
10.0 dB is 6.92 bits/s/Hz. 


Lastly, consider the 2x 2 matrix H given by 
el? -_e 
H= Se a } (2.88) 


This represents a MIMO system with N, = N, =2. The channel coefficients are fixed. 


The matrix Q is defined by H’’H and given by 
e og iO Toit _o-J8 
Petree. cid || oi git 


_/2-0 
‘es, Bf 
Then Equation (2.80) becomes 


C=W log, c C £1) 
Oo, 

C, P 

a 2 log, [5 


The capacity of the channel using this matrix H gives the same result as the pre- 


(2.89) 


(2.90) 


vious example when N =2. The capacity in this case also increases linearly with the 
number of antennas. If N =2 then the normalized capacity for the same SNR, 10.0 dB, 
is 6.92 bits/s/Hz. If the Q matrix can be reduced to a scalar multiplied by the identity 
matrix, then capacity is linearly proportional to the number of transmit antennas, maxi- 


mizing the possible capacity of the system. 
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D. THE MULTITPATH FADING CHANNEL 

The worst channel for any type of wireless digital communication system is the 
multipath fading channel. Performance in this type of channel is much less than that of 
AWGN alone [9,10,11]. In this section, an overview of the multipath fading channel is 
presented. The MIMO system design will be tested in this channel. Consequently, a 


model of this channel is derived in baseband equivalent form for use in the simulation. 


Firstly, radio waves may travel from a transmit antenna to a receive antenna by 
several different paths. This is called multipath. Multipath propagation of the signal is 
caused by reflectors and scatters present in the physical channel [9,10,11]. The signal, 
traveling via several different paths, will be sensed at the receive antenna several differ- 
ent times with different amplitudes, phases and arrival times [9,11]. This results in a re- 
ceived signal that can vary dramatically in amplitude and phase. In the frequency do- 
main, the spectral components of the signal may be affected differently by the channel. 
Thus, the frequency response of the channel may not be flat over the bandwidth of the 
signal. The overall result of multipath propagation is time—spreading of the signal that 


leads to intersymbol interference [9,10]. 


The configuration of the reflectors and scatters in the physical channel may 
change over time as well. This can be caused by the reflectors and scatters moving, or 
more likely, the transmit or receive antennas moving. As a result, the multipath channel 
is time—varying. The motion between the transmitter and receive antennas will cause a 


Doppler shift in frequency proportional to the relative velocity between them [9,10,11]. 


The time—spreading and time—varying aspects of the multipath fading channel can 
be measured and statistically estimated [9,10,11]. A maximum for the Doppler shift can 
be determined as well. Knowledge of this information aids the engineer in making key 


design decisions for a communication system. 


The time—spreading nature of the channel is characterize by the coherence band- 
width. “The coherence bandwidth is a statistical measure over which the frequency re- 
sponse of the channel is considered flat” [11]. The spectral components of a signal that 
fit in the coherence bandwidth will be passed with generally equal gain and linear phase. 


If a signal spectrum fits within the coherence bandwidth, then it undergoes flat fading. 
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Thus, the channel is flat fading or frequency—nonselective. If the spectrum does not fit 
within the coherence bandwidth, then it suffers frequency selective fading. Thus, the 


channel is frequency-selective [11]. 


The time-varying nature of the channel is described by the Doppler spread and 
coherence time [9,10,11]. The coherence time is inversely proportional to the maximum 
doppler shift [9,10,11]. The coherence time is the statistical measure of the duration dur- 
ing which the channel is essentially time-invariant. Thus, the coherence time gives an 
indication of how long the channel remains constant. If the symbol duration of a signal is 
less than the coherence time, then it undergoes slow fading. Thus, the channel is called 
slowly fading. If the symbol duration is larger than the coherence time, then the signal 


suffers fast fading and the channel is called fast fading [9,10,11]. 


The simulations were run for a slow fading, frequency—nonselective channel. To 
simulate a slow fading channel, the model of the channel coefficients is not permitted to 
change within a symbol period. This corresponds to the assumptions made concerning 
the channel coefficients in Sections A and B. In order to simulate the frequency— 
nonselective or flat fading nature of the channel, the channel coefficients must be distrib- 
uted according to a probability distribution function. Flat fading is generally simulated 
using a Rayleigh distribution for the magnitude of the channel coefficients [9,10,11]. 
This is consistent with a channel where the fading is caused by randomly moving scatter- 
ers and reflectors or moving transmitter and receiver. In order to achieve the baseband 
simulation of a Rayleigh distribution, two zero-mean, independent Gaussian random vari- 


ables are summed [11]. 


Let Z be the complex sum of two IID, zero-mean Gaussian random variables 


given by 
Z=X+jY =He*. (2.91) 
The mean and variance of X and Y are 
X=Y=0, (2.92) 


o_= Oo, =, (2.93) 


Let W be 


2 





W=|Z| , (2.94) 
or, equivalently 
See dae (2.95) 
The random variable W has a chi-squared distribution given by [10] 
ty (w) = | __y ferapwiret,, (w) , (2.96) 
o" Qn? T n 
2 
where the unit step function is denoted by 
0 w<0 
u(w)= 297 
( ) (i w=0. ( ) 


The variable n is defined as the number of independent random variables in the sum. In 


this case n=2. The distribution then becomes 


fy (w) = tenga yyy ; (2.98) 


This reduces to 
Ae (w) = oar (1) weer y (w) : (2.99) 


Since 
r(ij=(-)!=1, (2.100) 


the probability distribution function of W is given by 


Sy (w) =e"? u(w). (2.101) 
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The magnitude H of Z in Equation (2.91) is of interest and is given by 
=|Z|=VW. (2.102) 


The probability density function for W can be transformed to yield the probability den- 
sity function of H according to [12] 


Siu (h) = aa fair PO w) F. (2.103) 
where 
ee, (2.104) 
de 2(w) 


Using Equation (2.101) and the substitution defined in Equation (2.104) in Equation 
(2.103), we get 


a (h) = ahy? sre u(t) 


h (2.105) 
fx (A)= se u(h). 
o 
The mean and variance of H are given by [11] 
h=o a 
(2.106) 


Thus, H follows a Rayleigh distribution as required, and the sum of two zero-mean, ITD 


Gaussian noise sources will produce the desired channel model. 
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The phase © of Z is given by 











tan” Ea —1 for -o<x<0,-w< y<0 
xX 
ar 
© =, tan — for0<x<0,-w< y<oa (2.107) 
at x 
tan — +7 for -o<x<0,0< ys, 
Let A be 
Y 
A=—. 2.108 
7 (2.108) 


For the first case, when —0 < x < 0,-0 < y<0, the cumulative distribution function is 


given by [5] 
F,(a)=[" [fr (es y)dvae. (2.109) 


The probability distribution function is given by taking the derivative of F, (a) 





Hol 400 PY pa(xa)brde} oan 


By the Leibniz rule, the probability density function is given by [5] 


f,(a)=] 


0 


oS uy (x,ax) de. (2.111) 


The random variables X and Y are independent; therefore their joint probability density 


function is given by 
fr (9) = Fr OS (y). (2.112) 


The substitutions of Equation (2.112) as well as the Gaussian probability density function 
into (2.111) gives 





=) (on ; (2.113) 





Evaluating Equation (2.113), we obtain 











fala)= ae Yetta 
The random variable © is given by 
©=tan'(A)-z, 
where 
do|_ 1 
da| 1+a°’ 
and 


a= tan (0 + 1) : 
Using the transform defined in (2.103), we get 


1 


fo (@)=1+a On(l+a’) 





a=tan(0) 


Equation (2.118) simplifies to 


f,(0)=—- for “RSOS——. 


(2.114) 


(2.115) 


(2.116) 


(2.117) 


(2.118) 


(2.119) 


For the second case, when 0 < x <,-—0 < y<o, the cumulative distribution function is 


given by [5] 


rey (a) =| fe (x, y)dydx. 


(2.120) 


Like the first case, the probability distribution function is gained by taking the derivative 


of F,(a) 


dF, (a) 
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a STE Se (x.y)avde 


(2.121) 


By the Leibniz rule, the probability density function is [5] 
fle)=( Xf (x,ax) dx. (2.122) 


Similar to the first case first case, f,(a) is 














a)= o<a<on, 2.123 
f(a) 2n(1+a°) ae 
In this case © is given by 
@=tan'(A), (2.124) 
where 
do 1 
OG St 2.125 
da| 1+a° ( ) 
and 
a=tan(6). (2.126) 
After transformation f,(@) becomes 
1 a a 
0)=— for -—<O0<—. 2A27 
fo (9) Be 8 2 vee) 


For the last case, when —o0 < x <0,-0 < y <0, the cumulative distribution function is 


given by [5] 
0 a) 
F,(a)=[_ le fy (xy) dvds. (2.128) 
The probability distribution function is given by taking the derivative of F, (a) 


Hal) ATE Ff (asy)avar] me 





Once again, by Leibniz’s rule, the probability density function is [5] 


ioe) 


Sis (a)=J, —x fy (x, ax) dx, (2.130) 
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and f(a) becomes 














fal) Tse) o<sas<0. (2.131) 
In this case © is given by 
@=tan"(A)+z, (2.132) 
where 
< si — (2.133) 
and 
a =tan(0). (2.134) 
Finally, after transformation f,(@) becomes 
fo(0)=—— for =<0<x. (2.135) 
Dit Z 


Therefore, the complex sum of two zero-mean, HD Gaussian random variables 


gives a magnitude that is Rayleigh distributed with a phase that is uniformly distributed 


in [-x, x] . This models the baseband equivalent slow fading, frequency—nonselective 


channel. It will be used to simulate a multipath fading channel. 


E. DIVERSTIY AND ERROR CONTROL CODING 

As stated, the multipath fading channel is a severe environment in which to com- 
municate. The error rate performance suffers greatly compared to an AWGN channel 
[9,10,11]. However, two techniques available to mitigate the effects of multipath fading 
on performance are diversity and error control coding [11]. The possibility exists that a 
MIMO system can use both of these techniques in order to improve performance [1,13]. 
In this section, diversity techniques will be investigated. The Alamouti scheme, the focus 
of the simulations, incorporates diversity techniques in its space-time code, so the per- 


formance of diversity techniques is of great interest [1,13]. 
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The concept of diversity is to either transmit and/or receive the same information 
more than once. If one transmission suffers from a deep fade, then it is possible that 
other independent faded transmissions will not. The likelihood that all transmissions will 
be faded below possible recovery is reduced. The receiver, then, can choose from the 
best receptions [9,10]. There are three methods of diversity: time, frequency and space. 
Time diversity transmits the same information several times [9,10]. Frequency diversity 
transmits the same information on several carriers [9,10]. Lastly, space diversity em- 
ploys several antennas at the receiver in order to receive the same information several 
times [9,10]. In a MIMO system, multiple antennas are employed so space diversity is 
possible [1,13]. Also, frequency or time diversity is possible [1,13]. It will be seen that 
the Alamouti scheme uses both time and space diversity [1,13]. Therefore, the bit error 
performance of time and space diversity in a multipath fading channel will be considered 


for BPSK. 


1. Maximal-Ratio Combining 

The receiver that achieves the best performance in a multipath fading channel is 
one that uses maximal-ratio combining (MRC). The MRC receiver uses channel estima- 
tion, 1.e., CSI, of the magnitude and phase of each diversity reception. The MRC receiver 
multiplies the received diversity reception with the complex conjugate of the estimation. 
As a consequence, the phase is corrected for coherent detection and the magnitude of the 
signal is weighted by the signal strength. Strongly received signals will be weighted 
more than weaker signals (deep faded signals) in the decision process [9,10]. 

The time diversity reception of one of / =1,2,3,...L diversity signals for a BPSK 
system is illustrated in Figure 6. Each received diversity signal ,(t) = y,(¢)+n, (t) has 
duration 7, seconds. The correlation receiver equivalent of the matched filter is depicted. 
Assume that the signal after the correlation receiver is sampled at T, second intervals. It 
is then multiplied with the complex conjugate of the CSI represented by h,e™’” associ- 


ated with that time interval. 
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Figure 6. | BPSK with time diversity (After Reference [9].). 


Let x,, where k =0,1, be the transmitted baseband BPSK signal defined as 


YH xt) Aes SA for 0<¢<T, 
(2.136) 
x, =x,(t)=Ae” =-A for O0<t<T.. 
The energy F of a transmitted time diversity signal is [10] 
[, 2 T. 2 2 
E=( "(Of dt=J Of ae= 4°. (2.137) 


As the transmitted signal passes through the multipath fading channel in one time interval 


T. seconds, it is multiplied by a channel coefficient h,e’” and AWGN n,(t) is added. 
Thus, the received diversity signal y,(t) in one time interval is given by 
y(t) =hex, for0<t<T, , (2.138) 
where / is the time diversity reception and k =0,1. 
The symbol z;, represents the signal after the correlation receiver. The signal is 
given by 


2 = A he!"'x, +n,(t)dt 
AT he +N, when k =0 (2.139) 
eae +N, when k =1, 
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where N, is 
N,= Al n,(t)dt. (2.140) 


Substituting the energy E defined in Equation (2.137) into Equation (2.139), we get 


| |Ehe’” +N, when k = 0 
2 = (2.141) 
—Eh,e”* +N, when k =1. 
The signal z, after multiplication with CSI of the /-th reception is given by 
Eh; +he’"N, when k =0 
ees : (2.142) 
—Eh, +he’"N, when k = 1. 


The symbol z, represents the signal combined with all Z receptions. The signal z, is 
given by 
L . 
>| Bh; + he N, | when k =0 
l= 


ee (2.143) 
L 
-)>°| Eh +he’® N, | when k =1. 


I=1 
The comparator compares the real part of z, to zero. If Re|z,]>0, a binary 1 is 


the output. If Re[z, | <0, a binary 0 is the output. Then Re[z, | is given by 


> (en! +Re| h, e*'N, ]) when k = 0 
Rez 2s) (2.144) 
->'( Eh; +Re| he ven) when k =1. 
/=1 


The signal Re[z, | is a Gaussian random variable for a fixed set of {he \ [10]. The 


probability of an error P, is given by [9,10] 


ho) Zs } (2.145) 
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where Z* is the positive mean of the random variable Re[z, |, i-c., the expected value of 
Re[z, | when x, is transmitted and Orel is the variance of Re[z, |. The Q-function 


Q(x) is given by [5] 


ae 
O(x)=r= I, e dé. (2.146) 


Thus, to find P,, it is necessary to find Z* and Orel] The positive mean Z* is given by 


(2.147) 


The baseband equivalent noise n, (t) is complex sum of two IID Gaussian noise 


sources given by 
n,(t)=x,()+ jy,(0). (2.148) 


The autocorrelation functions of the baseband equivalent noise and its real and imaginary 


parts are defined as [5] 


R 
R,,(z) =E| x, (¢)x,(¢+7) | (2.149) 
R 


The power spectral density of the baseband equivalent noise, or its real and imaginary 


parts, is the Fourier transform of their respective autocorrelation functions. The two— 


sided baseband equivalent power spectral density S,, ( f ) for each diversity reception is 


defined as [10] 
B 
N, ls5 
S,, (f= 2 (2.150) 
0 fl> >> 
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where B is an arbitrary bandwidth that is B >1/T,. The power spectral density of the 


corresponding real bandpass noise S ( f ) centered about a carrier frequency f, is [10] 














N B 
g < 
r aes 5 
S(f)=4N B (2.151) 
0) pe 
ri If + ffs 5 
0 otherwise. 


The variance of the baseband equivalent noise O;, is equal to the total power of 


the noise given by [5] 
ao =N.B. (2.152) 
The total power of the baseband equivalent noise is the sum of the real and imaginary 


arts [5]. Let o2. represent the variance of the real part and o*. of the imaginary part. 
p nx, TEP p ny) ginary p 


Then the total power of the baseband equivalent noise is given by [5] 


O,, = 64, +o° (2.153) 


ny, * 


The real and imaginary parts of the noise are ID with equal variance. Then from Equa- 


tion (2.153) the variance of the real part Orn can be derived by 


o,, = 20%, 

207, =N,B (2.154) 
2 _NB 
nx) 2 


N 
a es 

S, (f)= . (2.155) 
0 fl>5- 


The noise power after the matched filter and multiplication by the channel parameter is 


given by [7,10] 
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0,7 =[ pier x,y af 








ae (2.156) 
-— [Lowy 4. 
Using Parsevel’s theorem in Equation (2.156), we get [5] 
2 
2 (2.157) 
_N AWE 
re 


The variance Orel is the sum of the noise power in each diversity reception and is given 


by [9] 

Pe tes x 2 
Re[z] = zy 

L 2 

=e 

NEX< 

el 
2 ‘A 








"led 
{ A (2.159) 


y=—Yh?. (2.160) 
Nae 
Then Equation (2.159) becomes 
P,=0(,27). (2.161) 


Note that 


_ Eh; 
V1 N, 





(2.162) 


represents the instantaneous SNR in the /-th reception [10]. 


In the course of developing the multipath model, the distribution for h; is already 


known. Ifthe random variable W, is defined as 
W,=H;, (2.163) 


then the probability distribution function is 


1 -w,/207 
fu, (,) =e 7'u(w,). (2.164) 


Note that the average SNR per diversity reception 7, is given by 








Eh? | E E E2o° 
7, =E|y,]=E}| —- |=—E| h? |= —E]w,]= 2.165 
Y1 [v7] = N, [ A N, [w, | N, ( ) 


The characteristic function of W, is given by [10] 


1 
F, (@)=————— (2.166) 
(2) (1- j207e) 
The probability distribution function of the sum of Z IID random variables is the 
L-fold convolution of the original probability distribution, or the multiplication of Z 


characteristic functions [5,10]. If W represents the sum of Z random variables W,, then 


the characteristic function of W is 


oe (2.167) 


(1- j20°w)" 


This is the characteristic function of the chi-square distribution of degree 2. Therefore, 


the probability distribution function is given by 
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L-1 


Ww _w/2e? 
Sy (w) “3eE(L)° u(w). (2.168) 


The parameter L will always be an integer so 

















I(L)=(L-1)!. (2.169) 
Then Equation (2.168) becomes 
we! is 
fA) = sre Emi’ wre"u(w). (2.170) 
Transforming this distribution with 
iS —_W F (2.171) 
0 
where 
E 
a =aee (2.172) 
we get 
1 uid /20° 
WMO-F are Tayi’ : (2.173) 
No wat 
7, 
Equation (2.173) becomes 
f(y)= = er (2.174) 
[f za (L-1)! 
0 


Substituting the average SNR 7,, Equation (2.165), into Equation (2.174), we obtain 


fe(y)=————e “u(y). (2.175) 


Equation (2.175) is the probability distribution function for T° as a function of 7. 
The probability of bit error is also written as a function of vy. The average probability of 


bit error in a multipath fading channel can be derived by taking the expectation of the bit 


error probability, that is [9,10], 


P.=ELA(7)] 
a (2.176) 
=[ By) f-(r)ar. 
Substituting Equations (2.161) and (2.175) into (2.176), we get 


L-1 a 


R= |) O(V27) amy 


dy. (2.177) 


This has a solution given by [10] 


7-4 a] si rs 5 4). (2.178) 


where the parameter wu is given by 





ae ee (2.179) 
1+7 


In time—diversity, the combiner adds L diversity receptions in succession. The 
transmit interval 7) is related to the bit duration 7, by [9] 


i irae BY Ae (2.180) 
Let the energy per bit Z,, be defined as 
EAT: (2.181) 


The correlation filter of the MRC integrates over one diversity reception 7,. The energy 


per bit relates to energy per reception by [9] 


E,=LE, (2.182) 
or 
E 
E=—, D183 
Fi ( ) 


Substituting Equation (2.183) into Equation (2.165), we get 





(2.184) 


In space—diversity, the MRC is configured as in Figure 7. The space—diversity re- 
ceptions come from separate receivers, and the time interval T of the matched filter is one 


bit period 7,,i.e., J, =7. Therefore, the energy per diversity reception is equal to the 


energy per bit 
E=E,. (2.185) 


The average SNR per reception becomes 





(2.186) 








aie 


A 


he 





Figure 7. | MRC space diversity configuration (After Reference [9].). 
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Of special interest is when L =1, i.e., no diversity. The average SNR per recep- 


tion for both time and space diversity schemes reduces to 





(2.187) 


P= | (2.188) 


Substituting Equation (2.179) into Equation (2.188), we obtain 


Pes te || (2.189) 
2; 1+y, 


This result is the probability of bit error P for a BPSK SISO system in a multipath fad- 





ing channel [9,10]. Equation (2.189) will be the benchmark for simulations in a multi- 
path fading environment. Furthermore, the means to analyze an MRC receiver is ob- 


tained. This will be essential in the analysis of the design. 


F. SPACE-TIME CODING 

“An effective and practical way to approaching the capacity of MIMO wireless 
channels is to employ space-time coding” [1]. Space-time coding is simply transmitting 
signals over the channel from multiple antennas to multiple antennas in such a way that 
capacity is maximized. Maximum capacity can only be approached. Space-time codes 
are designed to achieve both time and space diversity to improve the error performance 


and achieve diversity gain [1]. 


Another mechanism to improve performance is error control coding (ECC) [1]. 
Error control coding arranges the transmitted symbols in such a way as to increase im- 
munity to noise. It is realized by introducing redundant bits to the information bit stream, 
allowing the receiver to detect and possibly correct errors. However, increasing the 
number of bits in the bit stream increases the bandwidth required to transmit the informa- 
tion [10,11]. Ifthe error control coding is designed in conjunction with the modulation 
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technique and the transmit diversity scheme, then it is possible that there is no bandwidth 


expansion. With ECC, coding gain is achieved [1]. 


Overall, it is possible to achieve both diversity and coding gain with space-time 
codes [1]. (The general effects of diversity and coding on performance are illustrated in 
Figure 8 from Reference [13].). In general it depicts the possibility of improving per- 
formance by the introduction of diversity and ECC. The focus of this thesis is diversity 


gain; the addition of error control coding is left for future work. 


= Low SNR region 
~ = 


—_— 
Codinggain Diversity gain 


SNR (dB) 





SER 


Figure 8. _ Diversity and coding gain (From Reference [13].). 


In this chapter, the fundamentals necessary to understand a MIMO system in a 
multipath fading environment were introduced. It was seen that the capacity of a MIMO 
system can increase linearly with the number of transmit and receive antennas. To ap- 
proach this capacity, space-time codes are designed to incorporate diversity and error 
control techniques. A baseband model of a slow fading, frequency nonselective multi- 
path fading channel was obtained for simulation. Lastly, diversity techniques with the 
MRC receiver were examined to facilitate analysis of the design. In the next chapter, the 
Alamouti space-time coding scheme is introduced, and a transmitter and receiver design 


will be developed. 
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Hl. MULTIPLE-INPUT—MULTIPLE—OUTPUT TRANSMITTER 
AND RECEIVER USING THE ALAMOUTI SPACE-TIME CODING 
SCHEME 


This chapter introduces the Alamouti space-time code. Once the code is pre- 
sented, it will be used to design the MIMO transmitter and receiver using BPSK modula- 
tion. The designs will be simulated in Chapter IV and bit error performance will be ana- 


lyzed and tested. 


A, THE ALAMOUTI SPACE-TIME CODING SCHEME 
The Alamouti space-time scheme was introduced in Reference [14]. It is a 


space-time code for transmission through two antennas, NV, =2. It will be used ina 


MIMO configuration with two receive antennas, NV. =2. 


To describe signals more precisely from transmitter to receiver, the following no- 


tation will be adhered to through the remainder of this thesis. Let: 
e x” represents the signal x in the n-th transmission interval from the m-th 


antenna, 


m 


e x” represents the signal over the complete code sequence from the m -th 


antenna, 

e square brackets [ | encompass the total space-time code. The elements 
within the brackets separated by commas represent the code signal of one 
time interval, and 


e asterisk * represents the complex conjugate of the signal. 


Let x, and x, represent the baseband equivalent of signals for BPSK modulated 
bits over two consecutive transmission intervals of 7 seconds. Each is represented by 
+A. Let x' and x’ be the Alamouti space-time coded signals transmitted from antenna 
one and antenna two, respectively. Then the Alamouti space-time coded signals are de- 
fined by [1,13] 
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(3.1) 


Dee * 
x?=[x,,x1]. 


Although x, and x, are real in this case, the complex conjugate remains for generality 


and possible transformation to complex values in future work. 


The Alamouti space-time code, as it emanates from either antenna, contains in- 
formation necessary to demodulate successive information bits. The combination of the 
two consecutive signals from both transmit antennas corresponds to time diversity. Two 


antennas are used at the receiver, thus, space diversity is also achieved [1]. 


The Alamouti scheme has two important features [1]. First, the coded signals 
from the two transmit antennas are orthogonal over the code interval. If x' and x* are 


considered as complex vectors then they are orthogonal if their inner product equals zero 


[4,6]. Specifically, 
1(..2\'7 
xi(x’) =0. (3.2) 


Expansion of left hand side of Equation (3.2) gives 


x! (x?)" =| x, -*]/* 3.3) 
7 oe, - tx =), 


Thus, x' and x” are orthogonal. 


To describe the second feature, first let the code matrix be defined by X. The 


code matrix X is the combination of x' and x’ and is given by [1] 


x52 cs (3.4) 
x % x 
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Then, XX" is 


Ee x * * 
=i X. Xx. 
*L 1 2 2 
XX = . : 
Lx, x #5 x 
i * * * * 
XxX, + X4X, XX. — X5X, 
= * * * * 
| x5%, — XX, Ay Xy 1H Hy (3.5) 
[) 2 2 
Inf +pof 0 
7 2 2 
0 jal +h] 





(n° +ba) JIL. 


Equation (3.5) represents a scalar multiplied by an identity matrix. Note as well that the 
orthogonal feature of the Alamouti code is encompassed in this feature. Furthermore, 
this feature represents a condition similar to the example of the capacity of MIMO sys- 
tems in Chapter II. It was found that with this matrix, i.e., a scalar multiplied by the iden- 
tity matrix, the capacity increased linearly with the rank of H. This feature and its im- 


portance will become clear in the next chapter. 


With the Alamouti space-time coding scheme defined, the transmitter and re- 
ceiver can now be designed. A generic transmitter and receiver is described in Reference 
[1]. The transmitter and receiver described next is the implementation using BPSK 


modulation. 


B. TRANSMITER 

The objective of the transmitter is clear; produce two BPSK signals from an in- 
formation bit stream for transmission from two transmit antennas according to Equation 
(3.1). The transmitter design is shown in Figure 9. Let the binary information source be 
represented by 6. The binary information source undergoes serial-to-parallel conversion 
to produce two binary bits 5, and b, in successive intervals. Next, the two binary sym- 
bols are modulated by baseband BPSK modulators to give the corresponding two base- 


band equivalent signals x, and x,. In turn, x,, x, and their complex conjugates are 


routed to two parallel-to-serial converters to produce x' and x’ as defined by Equation 


(3.1). 
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An important feature of the transmitter is that the transmission rate of the infor- 
mation source R, is equal to the transmission rate of each of the two transmitted signals 
x' and x*. Each bit period is T seconds and is equal to the original bit period. The se- 
rial-to-parallel and parallel-to-serial conversion combination maintains the transmission 
rate of the binary information source R,. This is another important feature of space-time 
coding. In MIMO systems, time and space diversity can be achieved without bandwidth 
expansion[1]. In this case, the Alamouti space-time code achieves both space and time 


diversity while maintaining the original transmission rate R, through the channel. 


BPSK Parallel 









Modulator : 
Binary to Serial 
: Serial to x? 
Information 
5 Parallel 
intori BPSK Parallel 
to Serial 











Modulator 


Figure 9. | MIMO Transmitter. 


C. RECEIVER 
The transmitter was relatively straightforward to design. The receiver is some- 


what more complicated. Let the signal available at each receive antenna be the repre- 
sented by r' and r’. In the absence of AWGN and multipath fading (i.e., all channel 


coefficients equal one), r' and r’ are given by 


1 1 2 
r=x'+x 
(3.6) 
r=x'+x’. 


Incorporation of the definitions for x' and x* from Equation (3.1) into Equation (3.6) 


gives 
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’ ; : (3.7) 
r= ey | eam e 

Over the code sequence the same information is available at each receive antenna. This 

is possible because of the assumption that all channel coefficients are equal to one. Ina 

fading environment this will not be the case. Nevertheless, there is time and space diver- 

sity capability, an advantage that should be exploited. The logic of the receiver design is 


to exploit this and maximize the diversity reception of x, and x,. 


Let the signals be isolated in one transmission interval 7. A serial-to-parallel 


converter performs this function. The received signals by transmission interval are 
tata] * * 
r =| ny |=[x, +x. -% |, (3.8) 


and 





prea ery ee $ 
rath ty | =| Hey, x, |. (3.9) 


If the complex conjugate is taken on 7, and 7; the received signals become 


hax +x, 
ae (3.10) 
(7; ) aN NG 
and 
i =X, +X, 
: (3.11) 


From a closer look at Equations (3.10) and (3.11), there is a total of four signals 
available. The signals in each time interval, distinguished with different subscripts, rep- 
resent time diversity signals. The signals received by each antenna, distinguished by dif- 


ferent superscripts, represent space diversity signals. By combining correctly, both x, 
and x, can be isolated. Let z," and z;’, where m=1,2, be the time diversity combina- 
tion on the m-th space diversity reception that isolates x, and x, respectively. Then the 


time diversity combination for x, is given by 
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ze mee K” + (7;" ) 
=X, +X, +X, -X, (3.12) 


=2x,. 
Similarly, the time diversity combination for x, is given by 
va = rn” -(7") 


=x,+x,-(%,-x,) (3.13) 


= 2x,. 
Both time diversity combinations are available on both space diversity receptions. 
By combining these correctly, x, and x, can be further emphasized. Let z, and z, be 


the space diversity combinations that result in x, and x,, respectively. Then, z, is given 


by 


= 2x, +2x, (3.14) 


= 4x, 
and z, is given by 


7 Ae, 
= 2x, +2x, (3.15) 


= 4x,. 
From the logic of Equations (3.12) through (3.15), it is possible to build a receiver 


to take advantage of both time and space diversity receptions and isolate x, and x,. 


Lastly, the signals z, and z, need only be demodulated to produce b, and b,, 
where the ‘*’ distinguishes the demodulated binary signal from the original signal. 
Lastly, b, and b, are sent to a parallel-to-serial converter to reproduce the demodulated 
information bit stream 5. The receiver design is illustrated in Figure 10. 


The MIMO system that features the Alamouti space-time code is now designed. 
The system incorporates transmit and space diversity techniques in order to combat the 


effects of the multipath fading channel. In the next chapter, the design is implemented in 
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Matlab with Simulink and is tested for performance with AWGN alone and AWGN ina 


multipath fading environment. 
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Figure 10. MIMO Receiver. 
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IV. ANALYSIS AND SIMULATION 


In this chapter the error performance analysis and simulation of the MIMO system 
with the Alamouti space-time code is presented. The system, as described in Chapter III, 
was simulated in Matlab using Simulink. The system is analyzed and simulated in a pro- 
gressive manner from an AWGN channel to the multipath fading channel. To facilitate a 
smooth transition from each progressive step, the changes to the original simulation setup 


are highlighted. The simulation of the transmitter and receiver are described first. 


For ease of presentation, the Matlab simulation blocks, as they are used in the 
simulations, are presented in Appendix B. Further information concerning the blocks and 
the parameters are found in Reference [15]. The parameters, as depicted in Appendix B, 


are maintained throughout all simulations. 

Each simulation is compared to a baseband equivalent SISO BPSK system. An 
important parameter for comparison is the average energy per bit E,. In a single channel 
BPSK system, £,, is equal if a binary bit zero or one is transmitted. Let the modulated 
signal for a binary zero be represented by x,(t). If the amplitude over a period of T 


seconds is represented by A, then E, is given by [10] 
L. 2 2 
E,=(_ |x (¢)[ d¢=4°T. (4.1) 


The energy per bit is referenced to the output of the transmitter, thus the transmitted 


power into the channel for the SISO system is defined by 


E, 


a (4.2) 


Poso = 


The symbol £, is defined as the energy per symbol and is used to represent the energy 


per symbol period of the MIMO system. This is done in order to distinguish between the 


energies of the SISO and MIMO systems. The mathematical representation of E, will be 


derived later in this chapter. It will also be referenced to the output of the transmitter, 


and the transmitted power into the channel for one MIMO antenna is defined by 


a3 


Puuvo == - (4.3) 


A. SIMULATION OF THE MIMO TRANSMITTER AND RECEIVER 


1. MIMO Transmitter Simulation 

The block diagram of the transmitter simulation is given in Figure 11. The trans- 
mitter information source is a Bernoulli binary generator block. It produces the binary 
information source 6. The output is a binary bit stream of 1s and Os according to three 
parameters. The first is the probability of zero. Since the ones and zeros are equally 
likely, this parameter is set to 0.5. The second parameter is the initial seed for the ran- 
dom number generator. This value is a prime number that must be different from all 
other random sources in the simulation. By ensuring different seeds, the statistical inde- 
pendence of all sources is simulated [15]. The last parameter is sample time, which is set 


at the value 7 [15]. 
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Figure 11. MIMO Transmitter Simulation. 


The next two blocks in the MIMO transmitter are the buffer and deinterlacer. The 
buffer takes the input bit sequence b and produces a 2-by-1 frame at a sample time of 


2T seconds. The deinterlacer splits the frame into two separate outputs, isolating b, and 
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b, for manipulation. Together, the buffer and deinterlacer provide the serial-to-parallel 
conversion of the bit stream. The buffer must wait for 27’ seconds before producing an 
output that was generated by the Bernoulli binary generator. This causes a delay of 
2T seconds. 

Now that 5, and 6b, are separated, each is sent to a baseband equivalent BPSK 
modulator to produce the modulated symbols x, and x,. The modulation of the binary 


data is accomplished according to Table 1. 





Input 5,, | Output x,, 





0 A 





1 —A 
Table 1. | Modulation of binary data. 














where the amplitude 4=1. The modulated signal x, is routed through a complex conju- 
gate block to produce x,. Similarly, x, is routed through a complex conjugate block and 
a negative gain block to produce —x,. The baseband modulated symbols x,, x,, x, and 
—x, are then routed to the input of two interlacer blocks to form the final outputs. The 


interlacer and unbuffer blocks provide the parallel-to-serial conversion. The output of 
each unbuffer block has a period of 7 seconds. The output identified in Figure 11 as 
“Outl” is the code sequence for antenna 1, x'. Likewise, the “Out2” is the code se- 
quence for antenna 2, x”. Each conforms to the required code sequence of the Alamouti 


space-time code given in Chapter III. 


2. MIMO Receiver Simulation 
The block diagram of the receiver simulation is shown in Figure 12. The inputs 
indicated in the figure as “In1” and “In2” represent the received signals from antennas 1 


and 2, respectively, that is r' and r’. The buffer and deinterlacer blocks are identical to 


SD 


the transmitter and serve as the serial-to-parallel conversion of the received signals r' 
and r’. However, it is important to note, as with the transmitter, the buffer blocks intro- 
duce another delay of 27. The total delay through the system is now 47’. The math 
function blocks implement the complex conjugate of their respective inputs. The sum- 


ming blocks are the time and space diversity combining described in Chapter II. 
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Figure 12. MIMO Receiver Simulation. 


The baseband equivalent demodulator blocks demodulate the incoming signals 


according to Table 2 [15]. 





Input z,, | Output 5,, 





Re[z,,]J20| 9 








Re [zis <0 l 














Table 2. | Demodulation to binary data. 
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The output of the two demodulator blocks are the binary bit streams b, and b, . The final 
two blocks, the interlacer and unbuffer, serve as the parallel-to-serial converter producing 


the demodulated binary information b. 


B. SIMULATION IN AWGN 
The MIMO system is first subjected to AWGN only. In this simulation, it is as- 


sumed there is no fading. This implies that the channel coefficients 4, are all constant 


and equal to one. A block diagram of the MIMO system in AWGN is illustrated in 
Figure 13. Both the MIMO transmitter and receiver are collapsed into subsystems for 
ease of presentation. They remain the same as those described in Figure 11 and Figure 
12, respectively. In this simulation, two gain blocks are introduced at the outputs of the 
transmitter. This is done to normalize the total power transmitted into the channel for 
comparison to a SISO BPSK system. The MIMO channel in AWGN consists of the out- 
put of the two gain blocks to the two inputs of the MIMO receiver. The two AWGN 


channel blocks are complex white Gaussian noise sources that add to the incoming sig- 


nals. They are represented by n' and n’. 
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Figure 13. Simulation in AWGN. 
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1. Simulation Parameters 

The factor K in the gain blocks is determined from the discussion of the MIMO 
model in Chapter II. The power emanating from both of the transmit antennas is equal. 
Let the power from any antenna be represented by P*”. Let the total power transmitted 
into the channel be represented by P. The power from any antenna is the total power di- 
vided by the number of transmit antennas, in this case N, =2. The power from any an- 


Ant 


tenna P*” is given by 


Pt" = tes ee 
2 


P 
— 4.4 
N, ov 
The energy E, in any symbol x,, x, x, or x, is equal whether a binary zero or 


one is transmitted. The distinction between the SISO and MIMO system is the MIMO 
system symbols include the multiplication of K by the gain blocks. Once again, let the 


modulated signal for a binary zero be represented by x, (t) . If the amplitude over a pe- 


riod of 7 seconds is represented by A, then E, is given by [10] 
3 2 {pop go 
E, =| |Kx,(t)[@t=| K?#'dt= KAT. (4.5) 
The power emitted by each antenna is equal to the symbol power given by [7] 
fa ee eye (4.6) 
T 
When Equations (4.4) and (4.6) are equated, the total signal power into the channel as a 


function of amplitude and the gain factor K is given as 


a 2 42 
—=K’°A 
> (4.7) 
P=2Kea:, 
The total power at the input of either of the AWGN blocks with no fading is also 
given by P. This is set to 1.0 Watt. The amplitude A is 1.0 Volt. Using this informa- 


tion in Equation (4.7), we obtain 
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(4.8) 


Thus, the gain factor K is set to / V2 in this and all subsequent simulations. The energy 


per symbol £, of the MIMO system is 





(4.9) 


The energy per symbol for x,, x,, x, or x, is given by Equation (4.9). It incorporates 
the gain factor K. Thus, if a bit zero is transmitted in the MIMO system, represented by 


x,(t), it is given by 
x (t)=—  forO0<r<T. (4.10) 


Comparing the energy per bit LE, of the SISO system, Equation (4.1), and the energy per 


symbol E£, of the MIMO system, Equation (4.9), we have 
E, 

Eo=— 4.11 
5 (4.11) 


The AWGN channel blocks require several parameters: an initial seed, the energy 
per symbol-to-noise ratio E,/N, , the power of the input signal, and the symbol period. 
The symbol period is set to T seconds. The average energy per symbol-to—noise ratio 
E,/N,, increases with each simulation run. The initial seed is a prime number that is dif- 


ferent in each block as well as different from the initial seed of the Bernoulli binary gen- 


erator block. 


Finally, the error rate calculation block compares each bit from the original bit 


stream b with the demodulated bit stream 6 over the simulation run. It calculates the bit 


error rate or simulated probability of bit error. The total delay of 47 is incorporated in 


this block to ensure that the comparison of b and b is synchronized in time. 
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2. Performance Analysis 


: F - 2 . 
At each receive antenna, the received signals r' and r’ are given by 


(4.12) 


The complex valued vectors n' and n’ in Equation (4.12) represent the noise over two 
consecutive bits added at receive antennas | and 2, respectively. If the definitions for x' 


and x° from Chapter III are incorporated, then Equation (4.12) becomes 


r’ =([x.-%5 ]+[ 2524 |)+[ mm | 
r? =([,-x3 |+[ 42.27 |)+[ 42,72 |. 


After serial-to-parallel conversion, the signals in each time and space diversity reception 





(4.13) 





are given by 


s 

II 
— 
=e 
OL 
Li 

Il 


| (x, +2) +41,(41-25) +m |, (4.14) 
and 
r=[ee |=|(+e,)+m, (4-3) +m |. (4.15) 
The time diversity combination of the received signal on antenna 1, represented 
by z, and z,, respectively, are 
z=A+(a) 


* 


= (x +x,) +0) +((x7 - x3) +25) (4.16) 


= 1 * 
=2x,+n, +n, , 


and 


=(x +x,)+m -((x) -3) +4) (4.17) 


Similarly, the time diversity combination of the received signal on antenna 2, z, and z;, 


are 


z =2x,+m +n> (4.18) 


pia 2 2* 
2, SINAN, Ny 


The space diversity combining of x, and x, are represented by z, and z,. After 
space diversity combining, z, is 
aw Hee ae 
=2x, +n, +n, +2x,+n, +> (4.19) 
=4x,+m +n, +n +n; , 
and z, is 
Z,=2Z,+2Z, 
=2x,+nj—n, +2x, +n —n> (4.20) 
=4x,+n|—n, tn —n>. 
It is important to note the symmetry of the receiver. In Equations (4.19) and 
(4.20), x, and x, are the BPSK baseband modulation of two binary random variables that 
are IID. The noise components n,n, ,n,; and n; represent zero-mean IID complex 


Gaussian random variables. The sum of the means is zero, and the sum of variances is 
equal in both Equations (4.19) and (4.20) [5]. Thus, the probability of error in demodu- 


lating z, and z, are equal. The overall probability of bit error of the system is the aver- 


age of the two. Since they are equal, the probability of bit error of the system is the prob- 


ability of error in demodulating either z, or z,. 


Let the overall probability of bit error be represented by P,. To determine P,, z, 


is used. The demodulator takes the real part of z,. Let z be the real part of z,, given by 


z= Re(z,) = Re(4x, +n +n, +n; +n;') 
(4.21) 
= 4x, +Re(n|)+Re(n; )+Re(n;)—Re(n;’). 
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Let the sum of the real parts of all the noise components be represented by n. Then 


Equation (4.21) reduces to 
z=4x,+n. (4.22) 


Equation (4.22) defines a Gaussian random variable Z. The probability of bit error for a 


Gaussian random variable P, [9,10] is 





Z -}, (4.23) 
(oy 


Zz 


nad 


where Z* is the positive mean of Z (i.e., the mean of Z when a binary 0 is transmitted). 


Let %, (t) represent the modulated signal x, when a binary 0 is transmitted. If matched 


filter conditions are assumed, then the positive mean Z* at the output of the receiver is 


given by [9,10] 


=A" x3 (v)at (4.24) 


The power spectral density of the real part of each noise component is N, /2 
Watts per Hertz [15]. Let the total noise power of the real part of the noise be repre- 


sented by o?. The noise power co for each noise component after the matched filter is 





ae Vee N, 
Gal, 


a XA of. (4.25) 





By Parsevel’s theorem [5], Equation (4.25) converts to 


o, - iB - 





x(t) at 





(4.26) 


Let the total noise power of Z be represented by o?. Then o? is sum of each of the 


noise components co and is given by 
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a =40 SNE, (4.27) 


Using Equations (4.24) and (4.27) in Equation (4.23), we get 


me | AE. 
NE, 
: (4.28) 








Substituting the energy per bit defined in Equation (4.11) into Equation (4.28), we obtain 


of 


oO 


- (4.29) 
-o2 24 


oO 


The probability of bit error FP, for a baseband equivalent SISO BPSK system is 
given by [9,10] 








no? } (4.30) 
oO, 
In this case Z* is given by 
le (4.31) 
and 
2 = NA | (4.32) 
2 
Thus, P, is 


iy -o| |} (4.33) 


A comparison of Equations (4.29) and (4.33) reveals an improvement in bit error 
performance in the MIMO system by a factor of two over the SISO system. For equal 


performance, this equates to a 3-dB gain. 
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3. Performance Simulation Results 

The simulation was conducted with increasing E,/N,. The period T was set to 
0.0001 seconds. The simulation run time was set at 500.0 seconds. Thus, the total num- 
ber of bits in a single simulation run was five million bits. 

The simulated probability of bit error for the MIMO system is plotted in Figure 
14. The results are plotted with the theoretical results obtained in Equation (4.29). In 


addition, the theoretical probability of bit error for a baseband equivalent SISO system is 


plotted for comparison(Equation (4.33)). The simulated results follow the theoretical re- 


sults very well. 
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Figure 14. Results in AWGN. 

They deviate starting at approximately 5.0 dB and above, but only by approxi- 
mately 0.25 dB. This is a result of the number of bits in a single simulation run. If more 
bits are sent for comparison in a single simulation run, then comparison of simulated re- 
sults with the theoretical results improves. This was confirmed with several experiments. 
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Nevertheless, overall the simulation of the MIMO system with five million bits follows 


the theoretical analysis of bit error performance very well. 


C. SIMULATION IN A MULITPATH FADING CHANNEL 

Next, the effect of the multipath fading channels are added. The simulation block 
diagram is given in Figure 15. A multipath fading channel is inserted for each independ- 
ent path from transmitter to receive antenna. There are a total of four independent paths, 


each with a channel coefficient represented by h,,. 


1. Simulation Parameters 

The multipath fading model was derived in Chapter II. Each multipath fading 
model is identical with respect to configuration. A block diagram is given in Figure 16. 
The model is the complex sum of two IID Gaussian random noise sources with zero- 
mean and equal variance. To distinguish these Gaussian random noise sources from the 
AWGN, they will be called Gaussian channel fade noise sources. The three parameters 
of the Gaussian channel fade noise sources that must be determined and controlled are the 
variance, the sample time and the initial seed. The initial seed determines the independ- 
ence of the channel. Each of the seeds within the model is different. In this thesis the fo- 
cus is on the case of all channel coefficients being statistically independent. Therefore, 
each of the initial seeds of all the other models is different. Consequently, there are eight 


different seeds for the Gaussian channel fade noise sources. 


Next is the sample time. The sample time is determined such that the channel co- 
efficients remain constant over the code sequence. Thus, the channel is slowly fading 
over the code sequence. The duration of the code sequence is 27 seconds, thus the sam- 


ple time is set to 27 seconds. 
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Lastly, the variance of the Gaussian channel fade noise sources is set to a value 
such that the total average signal power entering both AWGN blocks is 1.0 Watt. To de- 
termine the variance of the noise source, the energy per symbol at the output of Multi- 


path! is examined. Let the energy per symbol at the output of the fading channel be rep- 


resented by E, given by [9,10] 


2 
_ tad 
B=) dt 





: (4.34) 
2 
aa] 47 
a5 
The channel coefficient is a random variable, so the average energy per symbol 


after the multipath channel is the mean of Z,. This can be obtained by taking the ex- 
pected value of Equation (4.34) [5]. Then E,, is 


8 a, 47 


. (4.35) 





2 
- FEL Ai 


The expected value of ln, is the second moment of the Rayleigh distribution. 


The second moment can be determined from the mean and variance of the Rayleigh dis- 


tribution defined in Chapter II by 
E| hy | =a) +h’ 
2 
2 (2-2 + 3 |Z (4.36) 
2 Z 
=20". 


Note that o* represents the variance of the Gaussian channel fade noise sources. Using 


Equation (4.36) in Equation (4.35), we obtain the average energy per symbol 


E,=o° AT. (4.37) 
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The average power per symbol P is 
PabaotA. (4.38) 


The signal power at the input to either of the two AWGN channel blocks is the sum 


power of two multipath channel blocks. Each of the channel blocks has equal average 
signal power. Then the total average signal power P. at either of the two AWGN blocks 
is 

P=2P=20°A’. (4.39) 


The total average signal power P must equal 1.0 Watt. The amplitude A is equal to 1 


Volt. Using this information in Equation (4.39), we find 


1=207 
- at (4.40) 
x 


Thus, the variances of all the Gaussian noises sources of the multipath fading channels 


are set to 1/2. 


2. Performance Analysis 
With the addition of the multipath fading channels, at each receive antenna the re- 
ceived signals are 
r= (hx! +h,x°)+n! 
: ; ; : (4.41) 
r = (yx +hy,x Jen 
Once again, incorporating the definitions for x’ and x” from Chapter III, Equation (4.41) 


becomes 


(4.42) 
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The channel coefficients do not change over the code sequence and the noise acts on one 
symbol period. After serial-to-parallel conversion, the signals in each space diversity re- 


ception and time interval are given by 














1 1 1 1 * i 1 
r= EF Fy =| Aim thx, +m Myx, — Nyx, + n; | : (4.43) 
and 
2 =| 7? |=|h h 2 hx, —hyx, +n 4.44 
POP oy P| Mar My Xy TM 5 My X — My Xq TF Ny | - (4.44) 
After time diversity combination of the received signal on antenna 1, z; is 
z=nt (4 ) 
T * ‘ay 
=h,x, +h x, +n, + (A, 2%, —h,,x,+ n;) (4.45) 
* * 1 1* 
= (hh, +h, )x, +(h,, =H, )%, +N +My, 
and z, is 


* 


1 1 1 
ay ire: (7; ) 
=x, +h,x, +n, —(Fo% Si, +n) (4.46) 
* * 1 1* 
=(M, —hiy) x +(h,; +h) Xx; +N — Ny. 
Similarly, the time diversity combination of the received signal on antenna 2 for x, and 


2 2 : 
x,, represented by z; and z; , respectively, are 


2) = (Ia, + Fig), + (Moy — A Jaa tmp tg” (4.47) 
22 =(hhy, — hy, ) x, + (Igy + hy, ) x, +0? — 
After space diversity combining, z, is 
2 = cA + Zz aa 


* * * 1 1% 2 2* 
= (Jy + iy + My + hy )X, + (Ay Ay + hy — hy ) xy +0 TN, tN +N; , 
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and z, is 


ed 2 
By 25 Zs 


, ‘ ‘ : ; . (4.49) 
= ( (Fy — ia + hay — Fy), + (Iya +A + Ig + yy) =) +] — tt, Ny 


Invoking the symmetry of receiver, we obtain the overall probability of error P. 
by using z,. The demodulator effectively takes the real part of z,. Once again, let z be 


the real part of z,. Then z is 


2=Re[z;| 
=(Re[M,]+Re[ ce ooetie hy, |)Re[] 
+(Re[h,,]—Re[ #, ]+Re[h,,]—Re[ h;, ])Re[x, ] 
ee ania al 


(4.50) 





The real parts of any h, are I[D, Gaussian random variables with zero-mean and 


equal variance. Their sum is also a Gaussian random variable with mean of zero and 
variance equal to the sum of the variances [5]. As already discussed, the real parts of the 
AWGN components are IID zero-mean with equal variance. Their sum is also a Gaus- 
sian random variable with zero-mean and variance equal to the sum of the variances. Let 


the sum of channel coefficients multiplied with x, be h,. Let the sum of the channel co- 
efficients multiplied with x, be h,. Let the sum of the noise components be represented 
by n. The real parts of x, and x, are equal to x, and x,. Equation (4.50) becomes 
zZ=h,x,+h,x,+n. (4.51) 
The overall probability of bit error, P,, can be determined by taking the sum of 


probability of error over all combinations of x, and x, [5]. If A’= A/ V2, then the real 


part of x, and x, over all combinations of b, and b, are 
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b, b, x 9) 
0 0 A’ A’ 
0 1 A -A' (4.52) 
1 0 —A’ A 
1 1 —-A' -A! 
Now P. is given by 
P, = Pr {error | b,b, = 00} Pr {b,b, = 00} 
+Prjerror|b,b, = 01} Pri bb, =01 
{error |b = 01} Pr{ybs = 01 — 


+ Pr {error | b,b, = 10} Pr {b,b, = 10} 
+ Pr {error | b,b, =11} Pr{5,b, = 113 


The probability b, and 6, being a bit zero or one are equally likely; therefore, 
Pr(b,b, = 00) = Pr(},b, = 01) = Pr(b,.b, = 10) = Pr(b,b, =11) = 1/4 and Equation (4.53) re- 


duces to 


Pie 7(Pr {error | b,b, = 00} + Pr {error | b,b, = 01} 


(4.54) 
+Pr {error | b,b, = 10} + Pr {error | b,b, = 11}) 
In terms of X, and X,, Equation (4.54) becomes 
P. = @r[z <0 X, =X, = A'|+Pr[Z <0|X, = 4',X, =-4'] 
4 (4.55) 


+Pr[Z >0| X, =-A',X, = A']+Pr[Z > 0| X, =X, =-4')). 
Concentrating on the first term of Equation (4.55) and incorporating this informa- 
tion in Equation (4.51), we get 
z=h A +h A +n. (4.56) 
Equation (4.56) represents a linear transformation of three independent Gaussian 


random variables h,,h, and n. The linear transformation of Gaussian random variables 


remains Gaussian [5]. Therefore, z is a Gaussian random variable with mean equal to 


the sum of the means and variance equal to the sum of the variances. Each of h,,h, and 


fal 


n have zero-means. Thus, z has a zero-mean. The first term of Equation (4.55) is given 


by [5] 


Pr[Z<0,X, =X, =4']=Pr[h,4’'+h,A'+n<0]=—. (4.57) 


N |e 


Each of the three remaining terms in Equation (4.55) gives the same result; thus, the 


overall probability of bit error is 


pe an ereites =e (4.58) 
4) 2 22> 2 2 





Therefore, the probability of bit error, P,, for the system is 0.5, regardless of the of the 


signal or the power of the noise. This was confirmed by simulation for increasing aver- 


age symbol energy to noise ratio E,/N, . 


3. Performance Simulation Results 


The simulation was conducted with increasing E,/N,. The period T was set to 


0.001 seconds. The simulation run time was set at 100.0 seconds. Thus, the total num- 


ber of bits in a single simulation run were 100,000. 


The simulated probability of error for the MIMO system is plotted in Figure 17. 
The results are plotted with the theoretical results obtained in Equation (4.58) In addi- 
tion, the theoretical probability of error for a baseband equivalent SISO BPSK system in 


a multipath channel derived in Chapter II is plotted for comparison. This is given by 


= Tt : Vy, 
R- |i fe) (4.59) 





where 7, is given by 





(4.60) 


Recall from Equation (4.40) that o” =1/2 thus, Equation (4.60) becomes 
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eas (4.61) 


The simulated results follow the theoretical results very well. However, they do not give 


the performance expected by a baseband equivalent SISO BPSK in a multipath channel. 
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Results in Multipath Fading Channel. 


The conclusion, of course, is that the receiver as designed in Chapter III and simu- 


lated in Figure 12 cannot communicate in a multipath fading channel. Thus, improve- 


ment to the receiver design must be incorporated in order for the system to function in a 


multipath channel. In the next section, the required improvement is incorporated in the 


design: the MRC receiver with CSI. 
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D. SIMULATION IN A MULTIPATH CHANNEL WITH MRC DESIGN 


An MRC receiver incorporates CSI [10]. The incorporation of CSI follows that 


found in Reference [13]. Nevertheless, before discussing the design of the receiver, the 


received signals are examined more closely. Investigation of Equations (4.43) and (4.44) 
once again gives 


1 
ry =h x, +hx, +n, 
* * 1 
r, =h,x, -—h,x,+n 
: hx, —h,x, z (4.62) 
Ko =hyx, +hy)x, +n, 


2 * ee 
Ty = Nyx, —hyx, +n). 


Next, we take the complex conjugate of 7; and r, , and Equation (4.62) becomes 


1 
a =h,x,+h,x, +n, 


1 
(7; ) = hyyx, — hx, +0, 
(4.63) 
te = hy x, +hyX, +1; 
2. 2* 
(7; ) = hyx,—hyx, +n). 


We would be well served to view Equation (4.63) as describing a 2-input, 4-output 


MIMO system. This equation set can be represented in matrix form as 


r=Hx+n, (4.64) 
where r is 
oi. Fi 
y 
r 
r= «|, (4.65) 
(72) 
(7) 








x 
x= } (4.66) 
xy 
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h, hy 
h h 
H= e 2 : (4.67) 
hy hy, 
h,, shy 
and n is 
ny 
n 
n=| "|, (4.68) 
ny 
2* 
LM 








Note that the second and third lines in Equation (4.63) were swapped in the definitions 
described in Equations (4.66), (4.67) and (4.68). 


From Chapter IT and Reference [1], the capacity is dependent on matrix Q given 


in Equation (2.74) by 


Q=H"H. (4.69) 


Using the definition of H in Equation (4.67), we get 


Q=H’"H 
h, hy 
-(h Toy. ti . hy tes (4.70) 
hy hy, Ni =H h hy, 
Vigo 
which reduces to 
(ls, °+[42f +f’ +P) 
oe Be cron aS (4.71) 
0 (lr, +|hg) +\-o,) + [to 


As discussed in Chapter II and Reference [1], this is a form of matrix Q that 
gives a linear increase in capacity with size of Q, in this case two. If CSI is used 


appropriately, then this capacity can be approached. In the receiver, if channel 
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priately, then this capacity can be approached. In the receiver, if channel information is 
used according to H” , then this can be achieved. 

It seems extraordinary that the matrix H has this property. However, it was the 
code and the reception of the code that forced this property [13]. Referring to Chapter 
III, we see that the code matrix X had this same property. Thus, in receiving this code, 


the channel matrix H has this property regardless of 4, , the channel realization over the 
code interval [13]. 
The use of CSI must be incorporated according to H’’. Given the received ma- 


trix r, let a signal matrix z be defined by 


* 


z=H"r. (4.72) 


This is consistent with Reference [13]. Thus, the implementation of CSI is accomplished 


according to Equation (4.72). 


The expansion of Equation (4.72) gives 





hy fy hy ty 








si his hy, h, —h, (": ) 
(1 ) (4.73) 
Z Aan thy +h, ie ) thy, (7 ) 
lis as hr, hy, (1 ) hay ( ) 
1. Simulation Parameters 


Equation (4.73) represents the rules for incorporating the CSI. The modified re- 
ceiver that incorporates the use of CSI is described in Figure 18. The modified receiver 
is a MRC receiver. The simulation of the MRC receiver is described in Figure 19. The 


channel parameters h,,,h,,,h,, and h,, are taken directly from the appropriate multipath 


fading block. This provides perfect channel knowledge. 
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MRC Receiver 


Multipath Simulation with MIMO MRC receiver. 
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In order to coincide with the correct symbols at the multipliers, the channel coefficients 
are delayed by 27 seconds since the buffers in the receiver induces this delay. The re- 


mainder of the simulation parameters remain the same as in previous simulations. 


2: Performance Analysis 


The received signal matrix z represents the signals z, and z, given by 


Zz * 
2-2 en, (4.74) 


Zy 
which expands to 


z=H’Hx+H'n 


: (4.75) 
z=Qx+H'‘n. 
After matrix multiplication, Equation (4.75) expands to 
(\a,. + |r) + lr.ol + lisa)’ ) x a hyn, a hyyny + hyn a hyn, 
-|* |- ate | 476) 
22 (\a, +|h,,| +|h,.| +|h,,| le —hyn, + hyn, —hyn; +hyn, 


Applying the symmetry of the receiver, we only have to investigate z,. The result 


of Equation (4.76) equates to an MRC receiver with ZL =4 diversity receptions [10]. The 


probability of error of z, was given in Chapter II, that is, 


ac ae 





where wu is 





w= |~L, (4.78) 
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and 7, is 





7; (4.79) 


The parameter E in Equation (4.79) represents the energy per symbol E,. The 


energy per symbol is given by 





E 
ye ae 4.80 
cs (4.80) 
Since o° =1/2,, Equation (4.79) becomes 
—. 
= . 4.81 
V1 2N, ( ) 
3. Performance Simulation Results 


The simulation was conducted with increasing E,/N,. As was the case in the 


AWGN simulation, the period, 7, was set to 0.0001 seconds. The simulation run time 
was set at 500.0 seconds giving the total number bits in a single simulation run of five 


million. 


The simulated probability of error for the MIMO systems is plotted in Figure 20. 
The results are plotted with the theoretical results given by Equations (4.77), (4.78) and 
(4.81). As well, the theoretical probability of error fora SISO BPSK system is plotted 
for comparison. Note that the SISO BPSK system has perfect knowledge of the CSI. 
The MIMO system’s simulated results, with MRC receiver, follow the theoretical results 
very well. The simulated results deviate starting at approximately 13.0 dB. This is a re- 
sult of the number of bits per simulation run. Overall, the simulation follows the theo- 
retical analysis of bit error performance. Also indicated in the figure is the improvement 


in bit error performance over the single channel BPSK system. For equal £,/N, , there 
is an excellent improvement in performance that increases as E,/N, increases. This im- 


provement follows the expected improvement as illustrated in Chapter II. Lastly, the sys- 


tem achieved L =4 diversity order; specifically, the sum of two time diversity and two 
79 


space diversity orders. This is consistent with the potential of the Alamouti space-time 


code [1,13]. 
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Figure 20. Results MRC Multipath fading. 


Found in Reference [1] are the results for probability of error fora MIMO system 
using the Alamouti scheme. They are given in Figure 21. The published results follow 
the same assumptions, i.e., the channels are independent and perfect CSI is known [1]. In 
the figure, the graph with the circle markings is the one of primary interest, which corre- 
sponds to the simulation in this thesis. The graph indicates very similar results, particu- 


larly at E,/N, =10 dB, where P. is approximately 10~. 


E. SUMMARY OF RESULTS 
In this chapter, the MIMO system designed in Chapter III was examined with 
AWGN only as well as the effects of multipath. The MIMO system performed as ex- 


80 


pected with only AWGN and achieved better performance than a SISO BPSK system. 
The system was than simulated for a multipath channel. As designed, it was not able to 
communicate through this channel. The receiver for the original system was converted to 
an MRC receiver with CSI. As a result, the new design was able to communicate suc- 
cessfully in the multipath channel. The performance of the new design matched expected 


theoretical performance. Additionally, it also outperformed a SISO BPSK system. 
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Figure 21. Published results (From Reference [1].). 


Furthermore, the system’s performance compared well with published results. The 
MIMO system was able to achieve L =4 diversity order with the Alamouti space-time 


code, which is consistent with published results. 


In the next chapter, the results of this thesis are summarized and areas for follow- 


on work are presented. 
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Vv. CONCLUSION 


The goal of this thesis was to develop a MIMO system with space-time coding 
with an eye towards military application. A MIMO system was designed using the 
Alamouti space-time code. It was simulated in Matlab with Simulink and was tested in a 
progressive manner, first in an AWGN channel and then in a multipath fading channel 
with AWGN. The performance of the system was analyzed and compared with simulated 


results, as well as published results found in Reference [1]. 


A. RESULTS 

The designed MIMO system performed in a manner consistent with theoretical 
analysis in the AWGN channel. The simulated results compared very well with the theo- 
retical performance analysis. However, the original design failed in a multipath channel 
with AWGN. The design did not incorporate CSI and a MRC receiver. Nevertheless, 
performance analysis was derived and simulated results reflected the anticipated theoreti- 
cal results. With this failure, the original receiver design was changed to an MRC re- 
ceiver incorporating CSI. Consequently, the new design was able to communicate suc- 
cessfully in the multipath channel. The performance of the new design matched expected 
theoretical performance, as well as published results. Finally, the designed MIMO sys- 
tem with the MRC receiver was able to achieve full diversity order with the Alamouti 
space-time code which was consistent with published results. There is improved bit error 
performance in the MIMO system over the SISO system and capacity gain which can be 


exploitable for military application. 


B. FOLLOW-ON WORK 

There are three areas identified for follow-on work. First, the designed MIMO 
receiver can be converted into a maximum-likelihood decoding receiver. This design 
holds promise of simplicity over the original design. The performance analysis can be 


investigated and simulated. 
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Next, higher order space-time codes can be investigated. Although the Alamouti 
space-time code can achieve full diversity gain, it is unable to achieve coding gain as 
discussed in Chapter II without bandwidth expansion. Higher order codes can be exam- 
ined and tested by augmenting the original design or using a maximum-likelihood decod- 
ing receiver. The incorporation of ECC can be investigated in conjunction with the 


modulation scheme in order to achieve coding gain. 


Lastly, the MIMO system is already using CSI at the receiver. If this information 
about the channel can be extracted at the receiver, then it can be sent to the transmitter. 
The transmitter can use CSI to more efficiently transmit power to the receiver through 
channels that are less faded than others. Techniques that control the transmitter can be 


investigated and performance of the system analyzed. 
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APPENDIX A. LIST OF ACRONYMS 





AWGN Additive White Gaussian Noise 
BPSK Binary Phase-Shift Keying 

CSI Channel State Information 

ECC Error Control Coding 

IID Independent Identically Distributed 
MIMO Multiple—Input—Multiple—Output 
MRC Maximal-ratio Combining 

SISO Single—Input—Single—Output 

SNR Signal-to-Noise Ratio 

SVD Singular Value Decomposition 
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APPENDIX B. MATLAB SIMULINK BLOCK PARAMETERS 


The block parameters shown in this appendix give the parameters and a brief de- 


scription of the Matlab blocks used in the simulations in this thesis. 














4 : wate 
[_ Frame-based outputs 


Samples per irame: 
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Block Parameters: BPSK Modulator Basebanal - 
BPSK Modulator Baseband [mask] [link] 
Modulate the input signal using the binary phase shift keying method. 


For sample-based input, the input must be a scalar. For frame-based input, 
the input must be 4 column vector. 


In case of frame-based input, the width of the output frame equals the 
product of the number of symbols and the Samples per symbol value. 


In case of sample-based input, the output sample time equals the symbol 
period divided by the Samples per symbol value. 


Parameters 
Phase offset [rad]: 


fo 
Samples per symbol: 


Block Parameters: Interlacer1 3 
Interlacer (mask) {link} 


Combine the elements of the input signals to generate the output signal. 
The elements of the first input signal become the odd-numbered elements 
of the output signal, while the elements of the second input signal become 
the even-numbered elements of the output signal. 


The inputs can be either scalars or frame-based column vectors. 


Block Parameters: Unbuffer1 


Unbuffer (mask) [link] 
Convert a frame to scalar samples output at a higher sample rate. | 


Parameters 
Initial conditions: 
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AWGN Channel (mask) [link] 
Add white Gaussian noise to the input signal. The input and output signals 


can be real or complex. This block supports multichannel input and output 
signals as well as frame-based processing. 


When using either of the variance modes with complex inputs, the 
variance values are equally divided among the real and imaginary 
components of the input signal. 


Parameters 

Initial seed: 

AWGN seed) 
Mode: [Signaltoncise ratio (Es/No) ss tttsti«iY 
Es/No [dB}: 
JE_No—“Osws—SSSSSCSSSS 
Input signal power [watts]: 
eff 


Symbol period [s): 


Block Parameters: BPSK Demodulator Baseband 
BPSK Demodulator Baseband [mask] [link] 


Demodulate the input signal using the binary phase shift keying method. 


For sample-based input, the input must be 4 scalar. For frame-based input, 
the input must be a column vector. 


In case of frame-based input, the width of the input frame represents the 
product of the number of symbols and the Samples per symbol value. 


In case of sample-based input, the sample time of the input is the symbol 
period divided by the Samples per symbol value. 


Parameters 
Phase offset (rad): 


jo 
Samples per symbol: 


———— 


OK | Cancel | Help | Apply | 
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Block Parameters: Gaussian Noise Generatorl 


ZMCG_seed(1] 








27 





Samples per Mame: 
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Block Parameters: Error Rate Calculation 


= 











| 


1_| 


Target number of errors; 








Mayimum MUMBEr oF symbols: 
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